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Abstract We discuss a consolidation of determinations of the density of neutral interstel- 
lar H at the nose of the termination shock carried out with the use of various data sets, 
techniques, and modeling approaches. In particular, we focus on the determination of this 
density based on observations of H pickup ions on Ulysses during its aphelion passage 
through the ecliptic plane. We discuss in greater detail a novel method of determination of 
the density from these measurements and review the results from its application to actual 
data. The H density at TS derived from this analysis is equal to 0.087 ±0.022 cm~ 3 , and 
when all relevant determinations are taken into account, the consolidated density is obtained 
at 0.09 ± 0.022 cm~ 3 . The density of H in CHISM based on literature values of filtration 
factor is then calculated at 0.16±0.04 cm" 3 . 

Keywords Heliosphere • Local Interstellar Medium • neutral interstellar gas • heliospheric 
interface • pickup ions • heliospheric UV glow 



1 Introduction 

With the plethora of detailed measurements available from heliospheric missions such as 
Voyager, Ulysses, SOHO, Cassini, ACE, Wind, and many others and in anticipation of the 
first mission dedicated to remote-sensing studies of the heliospheric interface IBEX, a con- 
solidation of heliospheric parameters is an important task. Key parameters in this quest are 
the density of neutral interstellar hydrogen and helium at the nose of the termination shock. 
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Whil e a consensus was reached on the parameters of interstellar helium already a few years 
ago dMobiuset al.[l2004l . 120051 : IWitteU 20041 : iGloeckler et all l2004al : IVallerga et all 120041 : 
lLallement et al.ll2004bl) . the hydrogen density remained more uncertain until very recently 
due to more complex effects on this species within the heliosphere. A dedicated team sup- 
ported by the International Space Science Institute from Bern, Switzerland, has recently 
reached a consolidation of the density of neutral interstellar hydrogen at the nose of the ter- 
mination shock based on different data sets and modeling. We discuss this consolidation, 
focusing mostly on the determination of the density at the termination shock based on hy- 
drogen pickup ion (H PUI) observations by Ulysses during its passage through the ecliptic 
plane at aphelion of its orbit, performed using a novel analysis method. 

2 Neutral gas density from pickup ion observations 

2. 1 Past attempts to derive the density at the termination shock 

Pickup ions are neutral interstellar atoms which become ionized and picked up by the frozen- 
in magnetic field of the solar wind. A simple model of the transport of pickup ions in the 
solar wind was developed by IVasyliunas and Siscod Jl976h under assumption that (i) the 
atoms of the seed (neutral) population are stationary with respect to the Sun at the moment 
of ionization; (2) the solar wind expands radially away from the Sun with a velocity Vsw, 
is spherically symmetric and stationary, and the frozen-in magnetic field is perpendicular 
to the expansion direction - thus there is no net parallel streaming of the PUI population 
(the PUI fluid expands radially); (3) the population of newly-injected pickup ions under- 
goes an immediate isotropization of the pitch angle with respect to the local magnetic field 
direction, forming a spherical shell distribution; (4) the evolution of the pickup ion distribu- 
tion function with distance is adiabatic; and (5) the ionization rate of the seed population is 
spherically symmetric and invariable with time. 

With these assumptions, the distribution function of pickup ions at a heliocentric dis- 
tance r at the upwind axis is described by the formula: 

where w = v/Vsw, v is the specific velocity of a pickup ion in the co-moving frame of the 
solar wind, and 

fton = Aon,E (r/r E ) 2 , (2) 

where jSi on .E is the production rate of pickup ions at the reference distance r E , usually 
adopted as equal to 1 AU. The production rate is assumed to be equal to the loss rate of 
the interstellar gas, which implies that tin is also a function of /3; on . 

It is evident from Eq.lQ~|i that the distribution function observed at a location (r, 6) for a 
given velocity v corresponds to the density of the seed population at a distance from the Sun 
equal to rw~ 3 / 2 . Thus, the absolute spectrum of pickup ions in velocity at a given location 
in space yields the density profile of the interstellar gas seed population with distance from 
the Sun if the PUI production rate and the speed of the solar wind are known. With the local 
density of the neutral gas in hand, one can extrapolate this quantity to the termination shock 
providing a sufficiently accurate model of the depletion of neutral interstellar gas in the inner 
heliosphere is available. Such a method has b een u sed to infer the densi t y of in terstellar H 
at the termination shock bv lGloecklej d 19961) and IGloeckler and Geissl J200ll) . A similar 
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method was used by iMobius et all dl988h to determine the termination shock density of 
interstellar He from observations at 1 AU. 

In some cases, an alternative method can be used. Since He + pickup ions are solely 
created by charge exchange between the neutral He atoms and solar wind alphas, and since 
Ulysses SWICS is able to simultaneously observe the He 2+ PUI a nd solar wind fluxes a t 
~ 5 AU, where neutral interstellar He is not significantly depleted, lGloeckleretal] ( ll997h 
could obtain the local He density with an uncertainty that depends only on the accuracy of 
the knowledge of charge exchange cross section and does not depend on the absolute cal- 
ibration of the instrument. This method, however, is not directly applicable to H because 
the production of pickup ions depends on two comparable ionization channels, i.e. charge 
exchange with solar wind and EUV ionization. Consequently, the rate cannot be easily elim- 
inated from the analysis by simultaneous measurements. The remaining uncertainty in the 
total ionization rate tran slates into the uncertainty in the density and its significance is il- 
lustrated by the fact that iGloeckler and Geissl d200lh had to invoke a high electron ioniza- 
tion rate of 2.4 x 10~ 7 s _1 at 1 AU to explain the pickup ion velocity distribution, while 
current estimat es for the time i nterval of the observations suggest that it is equal to about 
0.8 x lfr 7 s" 1 tezowsld.l2008h . 



2.2 Modeling uncertainties in pickup ion interpretations 



Because the aforementioned method to deduce the density at the termination shock depends 
on the knowledge of the neutral gas density inside the observer location and interplane- 
tary parameters greatly influence the relation between the local density and the density 
at the shock, it is rather sensitive to modeling choices and the knowledge of parameters, 
such as ionization rates and radiation pre ssure. The transport of PUI from t heir birthplace 
to the detector is qu ite complex (see, e.g. Ilsenberd. fl987l : IChalov and Fahj. 1 19981 : IChalovl , 
l2006l : lMalama et alll2006h and the attenuation processes of interstellar gas, including ion- 
ization and radiation pressure, ar e strongly aniso tropic and time-dependent on time scales 
from days to decades (see, e.g., IBzowskii [2008). Thus, for an instantaneous distribution 
of PUI (1) one cannot assume the local production rate of PUI is equal to the loss rate of 
the interstellar gas, which is one of the factors shaping «h (r, 0), (2) one cannot be sure 
where exactly the boundary of the cavity is located because it changes with time, and (3) 
one can use neither the cold model, nor even the classical hot model jFah lll978t[l979l : 
iThomasl [1971 IWu and Judge! 1 19791 : ITallement et all 1 19851) for the local density of inter - 
stellar gas. The reasons for this, partly discussed already bv lRuciriski and Bzowskil \ 19961) . 
are as follows. Firstly, H is subject to a c omplex behavior of the radiation pressure, which 
depends on the radial velocity of an atom l lTarnopolski and Bzowskll2008l : lBzow'skll2008h 
and is proportional to the net solar flux in the Lyman- a li ne. Its magnitude varie s with 
the phase of the solar cycle and over a solar rotation (e.g., IBzowskii GOOlbl 120081) . Sec- 
ondly, the io nization rate is composed of charge exchange with the solar wind protons, pho - 
toionization tezowskil2001bl : lAuchere et allboOSllBzowski et all 120081 : IBzowskii 120081) . 
and electron-impact ioniz ation fezowski et all l2008l ~ Bzowskil " 2008b - A ll three contribu- 
tors vary both with t i me |Bzowskill2001bl. 120081) a n d with heliolatitude jMcCqmas et all 



19991 : iBzowski et all 120031 : lOuemerais et all 120061) . IGloeckler and Geissl d200ll) have de 



termined the total ionization rate that includes all of these components from the slope of 
the observed PUI distribution. However, this method works under the stated assumption that 
at least over the observation period the PUI production rate is equal to the loss rate of the 
interstellar neutrals. Also it should be noted, the modeling of the PUI distributions for all 
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Ratios of H densities, various models to classical hot model 
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Fig. 1 Ratios of densities of neutral interstellar H computed using various models relative to the results of the 
classical hot model. Shown is the profile along the Ulysses-Sun line in the solar equatorial plane towards the 
Ulysses aphelion. Solid line: Moscow MC model (i.e., conditions at the termination shock non-Maxwellian 
and non-homogeneous); dotted line: radiation pressure sensitive to the radial velocity of the atoms (non-flat 
solar Lyman-a line shape adopted); dash-dot: full time-dependent treatment of the radiation pressure and the 
ionization rate. 



comparisons with observations, except if just the PUI flux at the observer location, i.e. at 
the PUI cut-off is used, is ti ed to the reasonable assump tion that the full adiabatic cooling 
in the expanding solar wind t lVasvliunas and Siscoelll976h is valid, and assumption that has 
not been independently tested with observations. Finally, there is a strong departure of the 
hydrogen distribution function at the termination shock from a shifted-Maxwellian due to 
the fi l tering processes in the he l iospheric interface region ( iBaranov et all[l998l:llzmodenovl . 
200 1 ; Izmoden ov et all l2005al ; iHeerikhuisen et al" , 2006; llzmodenov and Baranov. 2006). 

An illustration how strongly including or excluding different physical processes in the 
modeling can influence the simulated density profile is shown in Fig. Q] It shows simulated 
densities as a function of distance from the Sun in the equatorial plane towards the aphelion 
position of Ulysses obtained with three different modeling approaches normalized to the 
density profile from the classical hot model of the interstellar gas. 

As can be seen in Fig. [T] at heliospheric locations relevant to the Ulysses in-ecliptic 
measurements, various modeling strategies lead to differences in the predicted gas density, 
reaching a 30% level at ~ 5 AU and increasing towards the Sun. The magnitude of varia- 
tions shown here should be regarded as a lower boundary because the model used does not 
take into account the variations in the distribution function in th e source region, i.e. at the 
term ination shock, which w ere shown to amount to about 10% bv llzmodenov et all |2p05b) 
and llzmodenov et all d2008h . Therefore, the classical hot model can only be regarded as a 
qualitative representation of the true density, and the quality of the approximation is reduced 
with a decreasing distance from the Sun. The model distribution function derived from the 
classical hot model is also affected by the inherent deficiencies of this model, of course, and 
its direct use to the interpretation of PUI measurements is not recommended. 
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Ratios of H densities, effect of ji and j3 uncertainty 
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Fig. 2 Effects of changes in the absolute calibration of the ionization rate and solar Lyman-a flux on the 
neutral gas density profile. Shown are ratios of densities computed using the hot model for reduced fl and j8 
values over nominal values obtained from recent measurements. Solid line: p reduced by 25%; dots, dash-dot, 
dashes: /j reduced by 10%, 20%, and 30%, respectively; dash-dot-dot: both p. and j3 reduced, respectively, 
by 25% and 30%. 



Another factor relevant for the calculation of the absolute value of the H PUI distri- 
bution function is the accuracy of the heliospheric parameters such as the solar wind flux, 
solar EUV flux, solar Lyman-a flux, and the shape of the solar Lyman-alpha line. The ab- 
solute calibration of the solar wind flux measured by variou s instruments has been know n 
with an accuracy of ab out 25% (cf the O MNI-2 data base, iKing and Papitashvilil |2005), 
and an early review by iBzowskil feOOlah ). Similarly, the absolute calibration of the flux 
in the solar Lyman-alpha line has chang e d within the past 30 years by about 30% (see 
Vidal-Madiaj, il975l : IWoods and Rottmarl 1 19971 ; llbbiska et all l200d : IWoods et all [2000 ; 



Floyd et al 



2002) and even now is known to no better than 10%. The influence of uncer- 



tainties in these parameters on the absolute densities of neutral interstellar hydrogen inside 
the Ulysses orbit near ecliptic is shown in Fig. [2] 

As evident from Fig. [2] uncertainties in the absolute ionization rate of 25% translate 
into uncertainties of ~ 25% in the density of interstellar gas at Ulysses at ~ 5 AU, and an 
uncertainty of 30% in radiation pressure leads to a density uncertainty of 10% at 5 AU. 
Both uncertainties increase dramatically towards the Sun. All these uncertainties have to be 
included in the determination of the H density at the termination shock from PUIs in the 
inner heliosphere. 



2.3 Robust method to determine the density at the termination shock 

Recently, iBzowski et all d2008h have presented a novel method to determine the density of 
interstellar H at the nose of the termination shock using measurements performed at the 
boundary of the cavity of neutral interstellar gas in the inner heliosphere. 
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2.3.1 Simplified analytic derivation 



Based on the distribution function of H PUI as proposed by I Vasyliunas and Siscod Jl976h 
and repeated in Eq. (fl}, the phase space density at w = 1 (i.e., for locally injected ions) is 
directly proportional to the local source function of pickup ions S(r), defined as: 



S(r) = Aon "HO,Aon (>•))■ 



(3) 



Using a cold model of neutral interstellar hydrogen dFalu l 1 1 968l : lAxfordl . 1 19721) and assum- 
ing perfect compensation of the solar radiation pressure and gravity, the gas density n (r) on 
the upwind axis varies as: 



n (r) = no exp 
The local PUI production rate is equal to: 

S(r) = Aon wo exp 



VlSM 



Aon*? 1 



VlSM>" 



(4) 



(5) 



and scales linearly with the density at the termination shock. Now we seek a distance ro such 
that the source function does not depend on Aon* requiring that: 



_ =0 = «(r ,Aon('-o)) + Aon('-o) 575 

ofton "Aon 



(6) 



Combining Eq.© and Eq.© we get: 



With 



the distance ro becomes: 



n (r ) = -Aon (n>) 



dn 



<9w(ro,Ac 
d Aon 



d Ac 



'() 



Vjsm 



or r Q : 



Vism 



(7) 



(8) 



(9) 



in connection with Eq. {2). Inserting this result to Eq. 10, we find that the exact location 
where the source function is insensitive to the ionization rate is found precisely at the bound- 
ary of the heliospheric cavity, i.e. where the density is decreased to 1/e of the value at the 
termination shock. To calculate no from the locally observed PUI source strength S(ro) the 
ionization rate must be determined. However, as long as the observations are made close to 
ro moderate uncertainties in the knowledge of the ionization rate only result in rather small 
uncertainties in the determination of the terminatio n shock density np. In a comparison be- 
tween different simulations and Ulysses PUI data iBzowski et al.l 1 120081) showed that this 
method is also robust against uncertainties in the radiation pressure and against variations in 
the modeling approaches. 
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2.3.2 Simulations of the H PU I production rate at Ulysses 

Bzo wski et al. (2008) took a subset of the Ulysses H PUI data obtained over 13 months 
in 1997 and 1998 during the Ulysses passage through the ecliptic plane at aphelion and 
averaged the observed distribution function to eliminate fluctuations of the local heliospheric 
conditions. The resulting local production rate of H PUI at Ulysses was equal to 7.26 x 
10~ 10 cm" 3 s~ 3 . 

The simulations to caclulate the density of neutral hydrogen in the inner heliosphere 
were performed with the use of a combination of two models. First, the density of inter- 
stellar hydrogen at the nose of the te rmination shock was calculated using the Moscow MC 
model Baranov and Malama (1993) with solar wind alphas and interstellar He + ions in- 
cluded ( llzmodenov et all 120031) . This density was split into two populations, primary and 
secondary, and their respective bulk velocities and temperatures were registered. The pa- 
rameters of the two populations were subsequently handed over to the Warsaw kinetic time- 
and latitude-dependent model dRuciriski and BzowskiL 1 19951 : iBzowski et all Il997l 120021 : 
iTarnopolski and BzowskiL 120081) . which calculated a time series of densities and PUI pro- 
duction rates from the primary and secondary populations at 16 evenly-spaced moments 
during the observation interval at distances and offset angles coinciding with the Ulysses 
positions. The production rates were averaged over the appropriate time periods so that they 
cou ld be compared with the value obtained from observations. 

IBzowski et alJ d2008l) started with a test whether the handover of the calculations from 
the Moscow MC to the Warsaw kinetic model is robust. The point was that the distribution 
function of the neutral gas at the termination shock is not independent of the offset angle 
from the upwind direction, but the Warsaw model makes this assumption. This check was 
made by running both models with identical parameters, assuming that radiation pressure 
and ionization rate are averaged over the solar cycle, spherically symmetric, and invariable 
with time. The version of the MC model used in the simulations is stationary and axially 
symmetric. The Warsaw model was de-rated to similar assumptions and run with the pa- 
rameters of the heliospheric gas populations as predicted by the Moscow model at the nose 
of the termination shock. Finally, the densities predicted by both models along a radial line 
connecting Ulysses with the Sun were compared. It was discovered that for this geometry, 
discrepancies between the results begin at ~ 3 AU from the Sun and increase inward, but at 
the Ulysses locations, the two models agree within ~ 1%. Hence, using the static and axi- 
ally symmetric MC model to determine the parameters of the gas at the termination shock 
and handing these parameters to a time-dependent, full 3D Warsaw kinetic model was a 
rea sonable and unbiased approach. 

IBzowski et al.l ( 120081) studied in detail variations of the solar wind density and velocity 
as a function of time and heliolatitude as well as radiation pressure, photoionization and 
electron-impact rates for the observations time. They developed a time- and heliolatitude- 
dependent model of the charge exchange ionization rate and a spherically symmetric, time- 
dependent models of the photoionization and electron-ionization rates, as well as a time- and 
radial velocity-dependent model of radiation pressure, covering a large time span straddling 
the interval of Ulysses measurements, and implemented them all in the Warsaw code. With 
those in hand, they calculated the PUI production rates at Ulysses for the 16 observation 
periods and averaged the result. 

The robustness of this result was verified in a multitude of ways. First, tests of the 
robustness against modeling strategies were performed. The reaction of the simulated PUI 
production rate to a down-grading of the model from its full 3D time dependence was tested. 
First, the fully time-dependent ionization rate and radiation pressure were replaced with 
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Model/data ratios of H PUI production rates 
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Fig. 3 Model-to-data ratios of the source function of the H pickup ions at Ulysses as a function of reduc- 
tion factor of the net solar Lyman-alpha flux used in the simulation. The group of lines with filled symbols 
corresponds to the simulations performed with the density at the nose of the termination shock equal to 
0.095 cm~ 3 , the open symbols correspond to 0.1 17 cm -3 . Diamonds correspond to the nominal values of the 
ionization rate and radiation pressure dependent on v r , stars to the values of the ionization rate reduced by 
25% and radiation pressure dependent on v, , and squares and triangles to radiation pressure independent on 
v r with the ionization rate, respectively, on the nominal level and reduced by 25%. 



static values relevant for the 16 observation points. Then, the sensitivity of the radiation 
pressure to radial velocities of the atoms was switched off. Tests of the latitude dependence 
in the ionization rate were performed. 

With this result in hand, the gas density at the termination shock was computed from the 
relation: 

„ TS = « TS sim (10) 

where the index "sim" refers to a simulated quantity and "ob s" to an observed val ue. The 
magnitudes of the combined effects are shoewn in Fig. [3l(from lizowski et alj, l2008). which 
presents the combined variation of all modeling parameters in comparison with the obser- 
vations. It is evident that the calculated PUI production rate at Ulysses only weakly depends 
on these variations. 

Each of the simulations was repeated with departures of the ionization rate and radiation 
pressure from their nominal (measured) values. In particular, the radiation pressure was 
reduced by a factor of 0.9, 0.8, and 0.7 (i.e., by 10%, 20%, and 30%), and the ionization 
rate by 25%. The local PUI production rate changed just by a few percent in the case of the 
ionization rate variation. The effect of radiation pressure reduction was more pronounced, 
but still within ~ 15%. 

Finally, it was checked whether the conclusions hold when the parameters in the LIC 
are changed. The entire simulation program was repeated with a different set of parameters, 
which introduced the simulated density at the termination shock increased by ~ 25%; the 
most important difference, however, was a totally different breakdown between the primary 
and secondary populations. The simulated production rates increased almost exactly by the 
same percentage (as can be seen from the fact that the pairs of lines between the two groups 
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shown in Fig.[3]are almost perfectly parallel), proving that indeed, the PUI source function 
at Ulysses scales linearly with the density at the termination shock. 

In a final test, the simulations were repeated with a different coupling factor between 
the protons and neutr al atoms. It could be show n that replac ing the charge exchange cross- 
section formula from lMaher and Tinslevl l l 19771) with that bv lLindsav and Stebbingsl ([2005) 
does not visibly affect the results. 



3 Determination of hydrogen density at the termination shock 

To de termine the density of neutral H at the nose of the termination shock, iBzowski et alj 
d2008l) took the simulated values of the production rate at Ulysses relevant for the current 
calibration of the ionization rate and radiation pressure with their uncertainty bars (which 
correspond to the points for the reduction factor equal to 0.9 or 1.0 in Fig.[3]l and calculated 
wts from the observed PUI fluxes at the Ulysses aphelion. Since the modeling uncertainty 
was just a few percent, almost the entire error in the result is due to the uncertainty of 
determination of the PUI production rates from Ulysses observations (mostly the geometric 
factor of the instrument), which is known to about 20%. This analysis brought the density 
of interstellar H at the termination shock to 0.087 cm~ 3 ± 25%. 



4 Consolidation of the interstellar H density 

4. 1 Density of H at the termination shock 

Determinations of the neutral H density at the termination shock, obtained using several 
different observation methods, were recently compared and consolidated in the framework 
of an International Space Science Institute Team. These values were then used to infer the 
H density in the CHISM based on a global heliosphere modeling. 

As discussed in the previous section, a re-analysis of the Ulysses SWICS H PUI ob- 
servations resulted in a H density at the termination shock of «h,ts = 0.087 ± 0.022 cm~ 3 , 
where all observational, m odeling, and ex t ernal parameter uncertainties are included in the 
overall uncertainty of 25% lBzowskietaI] { |2008l) . The neutral H density at the TS was also 
determined using two independent methods and data sources. 

One of them was an analysis of the sol ar wind slowdown resulti ng from its mass-loading 
due to the accumulation of H pickup ions. iRichardson et alj J2008h took the solar wind ve- 
locity measured by Ulysses at ~ 5 AU from the Sun, when it was at the same heliolatitude 
as Voyager 2 at ~ 80 AU, and using a ID MHD model of the solar wind propagation, they 
calculated an expected solar wind velocity at Voyager assuming there are no pickup ions. 
The difference between the modeled and actually measured values, equal to 67 km/s, was 
adopted as the slowdown of the solar wind due to the "friction" of the solar wind against 
the inflowing neutral interstellar gas. Subsequently, they repeated the calculations of the 
solar wind speed at Voyager, but this time assuming three values of neutral interstellar hy- 
drogen density at the termination shock: 0.07, 0.9, and 0.11 cm -3 . Having compared the 
measured slowdown value of 67 km/s with modeling results, they arrived at a value of 
nH,TS — 0.09±0.01 cm~ 3 . With the observational uncertainty including potential variations 
in parameters, such as the ionization rate, and the modeling approach, the overall uncertainty 
for this derivation can be placed at ±0.02 cm~ 3 . 
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A third method was applied bv lPrvoretal](l2008h . who used UV backscattering obser- 
vations in an approach that does not require any absolute calibration of the observed UV 
intensities. They used data from Cassini located at 10 AU, and Voyager at ~ 90 AU from the 
Sun, obtained during the same time interval 2003 - 2004. Both spacecraft observed about 
monthly-periodic variations in the backscattering signal. These were attributed to the scat- 
tering of the Lyman-alpha radiation from active regions, which rotate with the Sun and thus 
illuminate the surrounding interstellar gas like a "lighthouse". The amplitudes of the varia- 
tions observed at Voyager were systematically lower than at Cassini, which was interpreted 
as attenuation by a factor of 0.2, resulting from multiple scattering of the Lyman-alpha pho- 
tons in the ambient neutral interstellar gas in the inner heliosp here. Using the Moscow MC 
model of the heliosphere and the radiation transport code by Gangopadhv av et all ([2006), 
they simulated the attenuation for two sets of interstellar parameters: (n p = 0.06 cm~ 3 , nn = 
0.18 cm~ 3 ) and (n p = 0.05 cm~ 3 ,nH = 0.15 cm -3 ), which brought the interstellar gas den- 
sity at the termination shock equal to 0.085 and 0.095 cm~ 3 , respectively. While it was 
impossible to exactly reproduce all observed variations, they concluded that the attenuation 
is consistent with a density at the termination shock of ~ 0.09 cm— 3. However, the value 
derived is model-dependent. Therefore, an overall u ncertainty of about 25% has been placed 
on the derived value of « h ,ts = 0.09 ± 0.024 cm" 3 dPrvor et al.ll2008h . 

Although aimed primarily at a determination of the ionizatio n rate in the inne r helio- 
sphere, another recent study of heliospheric backscatter glow by Ouemerai s et al ] ( l2006h . 
which is based on multi-year photometric all-sky observations with SOHO SWAN, has also 
returned a value of 0.09 cm~ 3 for the H density at the termination shock. 

Hence, all three results obtained using independent methods and independent data sets 
are consistent with each other with in their uncertainty bands, and with the result of the 
PUI analysis by Bzowski et al. (2008). Computing a weighted average from all these values 
as done for He by Mob ius et al.l ( 12004) . the density of neutral interstellar hydrogen at the 
termination shock is found equal to 0.089 ± 0.022 cm -3 . 

As an alternative attempt, Gloeckler et al. (2008, this volume) combined the CHISM 
density of neutral He, based on the c ombined He parame ter determination with a density 
"He = 0.015 ±0.0015 ±0.002 cm" 3 ( iMobius et all l2004h . with the abundance ratios ob- 
tained with Voyager from the suprathermal tails of interstellar pickup ion distributions and 
from anomalous cosmic rays in the heliosheath. They present a noticeably lower H density 
at the termination shock of 0.055 ±0.021 cm~ 3 , but with a rather large uncertainty due to a 
number empirical modeling steps involved in the method. 



4.2 Extrapolating the termination shock density to the CHISM 

The bulk velocity and temperatu re of the gas in the CHISM are now well established from 
interstellar helium obs ervations JWitteL 2004 ;lGloeckleretal.ll2004bl : IVallerga et all |2004 
Lallement et al ., 2004a), and compiled bv lMobius et al.l (2004), but the density, ionization, 
and magnetic field still require further studies. While the determinations of the density of in- 
terstellar H at the termination shock based on various observations, discussed earlier in this 
section, depend relatively weakly on the modeling of the heliosphere as a whole, extrapolat- 
ing this density to the CHISM is certainly model-dependent and requires a detailed model of 
the interaction of the solar wind with the local interstellar medium. Since the mean free path 
of the interstellar H atoms is comparable with the size of the heliosphere, kinetic models are 
required. In fact, due to modifications in the interaction region the flow of H atoms at the 
termination shock varies with the angle off the upwind direction. This dependence should be 
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taken into account in the analyses. In this paper we assumed that the hydrogen interstellar at 
the termination shock gas consists of two po pulations (primary an d secondary) which flow 
homogeneously at the TS. It was verified (see lBzowski et alll2008h that such an assumption 
does not introduce an error greater than 1% at the location of Ulysses observations discussed 
in this paper. 

Under these assumptions we c an estimate t he H atom number density in CHISM by in- 
troducing the filtration factor (e.g. Ilzmodenovl . 120071) as: F = «H,Ts/ n H,CHiSM, where «h,ts 
is equal to the density of neutral interstellar hydrogen at the nose of the termination shock 
and wh.chism to the H density in the unperturbed interstellar gas. Then, knowing «h,ts up- 
wind and the filtration factor from the models, we can obtain an estimate for mh.chism- The 
filtration factor should be calculated in the framework of a solar wind/CHISM interaction 
mo del. 

Ilzmodenov et al.l 11200411 present an axisymmetric stationary kinetic -gasdynamic model 
of the solar wind/CHISM interaction. The model differs from the original Baranov and Malama 
dl993h by taking into account the solar wind alpha particles and the ionized interstellar he- 
lium co mponent. The solar w ind param eters in this parametric s tudy were similar to those 
used by lBzowski et al.l J2008h . Table 1 in llzmodenov et all d2004h presents the filtration fac- 
tors obtained for different interstellar proton and H atom number densities. It is shown that 
the filtration factor varies from 0.51 to 0.58 for a wide range of densiti es. Applying t his fil - 
tration to the density at the nose of the termination shock obtained by lBzowski et al.l {2008) 
we obtain «h,chism = 0. 16 ± 0.04 cm~ 3 . However, since the error bars for wrts are large, a 
wider range of («h,chism, «p,chism) is needed than available from Izmo denov et al.1 d2004h . 

In addition, filtrat ion factors also depend o n the interstellar magnetic field, which was not 
take n into account byllzmodenov et al.l ( 120041) . a nd on variations in the sol a r wind ram pres- 
sure (McComas et al.l ( 120081) : for model, see, e.g. Jlzmodenov e t al. ( 2005ah, [izmodenov et al.l 
d2008n T Also, the dependence of the solar wind on heliolatitude may be important. Includ- 
ing these effects may extend the range of the filtration factor beyond the values mentioned 
above. To establish the interstellar parameters based on w H ts, additional constraints, such 
as the location of the TS in Voyager 1 and 2 directions (see, e.g., Izmodenov, this volume), 
Lyman-alpha absorption (Wood et al, this volume), and solar backscattered Lyman-alpha 
radiation observations (Quemerais et al., this volume), should be added. 

It is worthwhile to mention that the interaction of int erstellar gas with the solar wind is 
also described with hydrodynamic multi-fluid models (e.g.lMuller et all l2006: Sche rer and 
2003; iFlorinski et all I2005L 120031) . lAlexashov and Izmodenov ( 2005) compared results of 
their kinetic model with different multi-fluid models. Results relevant to present studies are 
listed in Table 2 of this paper and show that the parameters of the primary and secondary 
populations of interstellar atoms obtained from hydrodynamic models differ significantly 
fro m the results of the k inetic models. 

iMiiller et~aH d2008h compared the results from five different models of the heliosphere 
evaluated with identi cal boundary conditions. Two of these models are fully kinetic in their 
treatment of neutrals tearanov and Malamal 19931: Heerikhuisen et al. L 20061) . and three are 



hydrodyna mic mu lti-fluid dFlorinski et al.U2003[|2005l) : lMuller et al.U2006h . extended from 
IPauls et alj l ll995l) . andlScherer and Fahrl d2003h . extended from lFahr et alj | |2000|) ). The pa- 
rameter set used bv lMiiller et al.ld2008h was clos e to one of the sets of interstellar parameters 
used in the simulations by Bzowski et al.l d2008h with the important exception that the solar 
wind speed was essentially lower. The filtration factors obtained in this comparison varied 
fro m 0.52 (for the two -fluid hydrogen model) to 0.69-0.74 (for the kinetic models). 

Miill er etall d2008l) did not present the parameters of the primary and secondary popu- 
lations at the termination shock separately and offered only the mean values of the density, 
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bulk velocity, and temperature. While all models reproduced the mean bulk velocity at the 
termination shock within a narrow range from 19.2 to 23.4 km/s, differences in the temper- 
ature were larger, i.e. from 12 000 to 30 900 K. 

Since the goal of the comparison made bv lMuller et ai] ( 120081) was checking the results 
of different calculating schemes, all models had to include the same physical effects, which 
required switching off some of their advanced features (different in different models). As a 
result, none of them was ru n in its full mode. Nevertheless, both of the comparison studies 
dAlexashov and Izmodenovl d2005l) and iMiiller et al.1 d2008l) ) seem to suggest that an un- 
critical use of hydrodynamic multi-fluid models for quantitative analysis of heliospheric 
observations may be misleading. 

From the discussion above it follows that without performing a dedicated modeling, 
the H density in the CHISM can be assessed from the density at the termination shock by 
dividing this value by the filtration ratio being in the (0.51, 0.59) range, which brings the 
density of neutral H in the CHISM equal to 0.16±0.04 cm~ 3 . 

5 Conclusions 

The density of neutral interstellar gas at the termination shock seems to be well constrained 
by observations performed using different techniques and appropriate modeling. The datasets 
used include Ulysses observations of H pickup ions, Voyager and Ulysses observations 
of solar wind slowdown, Cassini and Voyager observations of modulation in the helio- 
spheric Lyman-alpha glow, and multi-year all-sky photometric observations of the glow 
by SWAN/SOHO. The densities obtained from analysis of these observations group in the 
range of 0.08 - 0.09 cm~ 3 and the uncertainty of this determination is about 20%. Extrap- 
olation of this value to the CHISM based on analysis of heliospheric filtration factors from 
the literature brings the density of neutral interstellar hydrogen equal to 0.16±0.04 cm~ 3 . 
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